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Objectives: We investigated the association between particulate matter less than 10 urn in aerodynamic diameter (PMio) exposure 
and non-accidental mortality in Asian populations by meta-analysis, using both time-series and case-crossover analysis. 
Methods: Among the 819 published studies searched from PubMed and EMBASE using keywords related to PM10 exposure and non- 
accidental mortality in Asian countries, 8 time-series and 4 case-crossover studies were selected for meta-analysis after exclusion by 
selection criteria. We obtained the relative risk (RR) and 95% confidence intervals (CI) of non-accidental mortality per 10 ug/m 3 in- 
crease of daily PM10 from each study. We used Q statistics to test the heterogeneity of the results among the different studies and 
evaluated for publication bias using Begg funnel plot and Egger test. 

Results: Testing for heterogeneity showed significance (p<0.001); thus, we applied a random-effects model. RR (95% CI) per 10 ug/m 3 
increase of daily PM10 for both the time-series and case-crossover studies combined, time-series studies relative risk only, and case- 
crossover studies only, were 1 .0047 (1 .0033 to 1 .0062), 1 .0057 (1 .0029 to 1 .0086), and 1 .0027 (1 .001 0 to 1 .0043), respectively. The non- 
significant Egger test suggested that this analysis was not likely to have a publication bias. 

Conclusions: We found a significant positive association between PMio exposure and non-accidental mortality among Asian popula- 
tions. Continued investigations are encouraged to contribute to the health impact assessment and public health management of air 
pollution in Asian countries. 
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INTRODUCTION 

Associations between air pollution and daily mortality out- 
comes have been well-documented by numerous studies [1- 
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4], and it has been reported that, among the ambient air pol- 
lutants, particulate matter (PM) has shown the most consis- 
tent associations with mortality [5-7]. Thus PMio (defined as 
PM less than 10 urn in aerodynamic diameter) concentration is 
monitored in most countries. As epidemiologic evidence is 
crucial for managing air pollution control-related policies, 
studies of the health effect of PM on the general and suscepti- 
ble population at large are ongoing. 

However, most air pollution studies and their recommenda- 
tions have been conducted in developed countries such as 
North America and Europe. Therefore, there is a need for stud- 
ies designed for Asian countries, with their different meteoro- 
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logical and socio-demographic conditions, as has been em- 
phasized in Health Effects Institute (HEI) report [8]. Much at- 
tention has been focused on fast-developing Asian cities with 
rapid industrialization, urbanization, and motorization, and at- 
tempts are being made to assess the health impact of air pol- 
lution by considering their unique characteristics. The HEI has 
also reported that PM may also be associated with total or 
non-accidental mortality in Asian cities of the more developed 
countries. 

Since its introduction in 1990, time-series analysis has been 
the method of choice in air pollution epidemiology. This meth- 
od examines the same population repeatedly under various 
exposure conditions, thereby allowing for time-invariant con- 
founders, such as age and diet, to be adjusted for. Also, the in- 
troduction of generalized additive models has enabled adjust- 
ment for the confounding effects of trends and seasonality 
and has boosted research on air pollution and health effect as- 
sociations [9]. However, time-series analysis has also been criti- 
cized for model complexity and the possibility of producing 
incorrect standard errors upon correlation between nonlinear 
functions [10,11]. 

Case-crossover analysis, designed to investigate the associa- 
tion between intermittent short-term exposure and rare acute- 
onset events within close temporal proximity, provides an al- 
ternative Poisson time-regression model for studying the 
short-term effects of air pollution. It has the advantage of con- 
trolling time-invariant individual confounders, and by match- 
ing case and control days within a short interval of usually 
within 28 days or a month, time-variant confounders are also 
removed. Thus, both time-series and case-crossover analyses 
are commonly used and are viewed as competing methods in 
air pollution epidemiology [1 2]. 

Thus, in this study, we attempt to investigate the association 
between PM10 exposure and non-accidental mortality in Asian 
countries by including both designs in searching the literature. 

METHODS 

Study Selection 

Literature searches were conducted using a combination of 
three groups of keywords: Asian country names (Korea, Japan, 
China, Mongolia, Thailand, Malaysia, Indonesia, Vietnam, Pap- 
ua New Guinea, Philippines, Singapore, and Taiwan); PMio, as 
exposure; and non-accidental mortality, as outcome. Thus, us- 
ing PubMed, EMBASE, and the Cochran Library, the searched 



keywords were [Korea OR Japan OR China OR Mongolia OR 
Thailand OR Malaysia OR Indonesia OR "Vietnam" OR "Papua 
New Guinea" OR Philippines OR Singapore OR Taiwan AND 
(PM10 OR "air pollution" OR particulate) AND mortality]. The 
study collection included all published papers until August 
201 2. Two authors (HYP and SB) independently evaluated the 
eligibility of all the studies searched by the keywords from the 
databases based on predetermined selection criteria, and 
agreement was found between the evaluators. In all the stud- 
ies, daily PMio was measured at public monitoring stations, 
and 24-h concentrations were provided. Non-accidental mor- 
tality data was obtained using the International Classification 
of Diseases (ICD)-9 (001-799) or ICD-10 (A00-R99) codes re- 
ported by a government or district institution. We only chose 
those results without co-pollutant effects. In the time-series 
studies, lag days of 0 to 1 (2 day average of lag 0 and lag of 1 
day) were selected a priori according to the literature [13,14]; 
if a study did not have results with a lag of 0 to 1 days, results 
with lag days of 0 to 2 were chosen as a second option, fol- 
lowed by lag of 1 day. 

Statistical Analysis 

We obtained the relative risk (RR) and the 95% confidence 
intervals (CI) of non-accidental mortality per 1 0 ug/m 3 increase 
of daily PMio from each study. Our main analysis included eight 
studies representing each study region, and we also performed 
separate meta-analyses with six time-series and four case- 
crossover studies. For sensitivity analysis, we performed meta- 
analyses with all ten of the chosen studies and compared the 
overall estimate with the main results in eight studies. 

We used Q statistics in a log odds scale with unequal weight 
for testing the heterogeneity of the results among the differ- 
ent studies, and considered a significance level of <0.1 as re- 
jecting the null hypothesis that the studies are homogeneous. 
We chose a fixed effects model on observing the homogenei- 
ty among the studies, and a random effects model on observ- 
ing the heterogeneity. Evaluation for publication bias was 
performed using Begg funnel plot and Egger test. For both 
tests for heterogeneity and publication bias, longitude and 
latitude information were used as modifiers. Meta-analyses 
were carried out using R software version. 2.14.2 "metaphor" 
package (Wolfgang Viechtbauer, Maastricht, The Nether- 
lands). 

Throughout the study process, we referred to the Meta-anal- 
ysis of Observational Studies in Epidemiology guidelines. 
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RESULTS 

Of the 81 9 studies searched by the keywords, 1 86 duplicates 
were removed, after which those studies not presented in 
English or not on human subjects were filtered. Studies with 
exposures other than PM10 (e.g., suspended particulate matter 
(SPM) and other air pollutants such as N0 2 , S0 2 , CO, 0 3 ) and 
mortality outcomes that did not include non-accidental mor- 
tality were excluded. 

The 48 filtered studies were reviewed thoroughly, including 
the entire text, and the study details were recorded. The title, 
authors and other details of publication, study period and city/ 
region, statistical method used, adjusted covariates, stated lag 
days, presentation of the results, geographical position (in lon- 
gitude and latitude), and population size, if available were in- 
cluded. Further exclusions were carried out using the follow- 
ing criteria: other than time-series or case-crossover analysis 
(n=8), unfit population (n=7), insufficient accuracy (n=3), or 
overlapping dataset, period, or region (n = 18). In the case of 
selecting studies from overlapping datasets, periods or re- 
gions, studies with results from more recent and longer peri- 
ods were selected. In total, eight time-series studies [14-21] 
and four case-crossover studies [22-25] were finally selected 
for the meta-analysis (Figure 1). 

Table 1 shows the main characteristics of the studies includ- 
ed in the meta-analysis for PMio-mortality association. We se- 
lected eight time-series studies (Wuhan, China; Shanghai, Chi- 
na; Seoul, Korea; Incheon, Korea; Hong Kong, China; Bangkok, 
Thailand; Guangzhou, China; andTianjin, China) and four case- 
crossover studies (Shanghai, China; Anshan, China; Seoul, Ko- 
rea; and Kaohsiung, Taiwan) for separate group analysis. For 
estimating the combined effect of both time-series and case- 
crossover studies, we excluded studies by Bae and Park [17] 
and Kan and Chen [22] as there was an overlap of region 
(Seoul and Shanghai, respectively). Studies by Yi et al. [24] and 
Kan et al. [16] were selected over the two excluded studies 
due to prioritization for recent data analysis. 

Before the main analysis, we tested the results from the 
studies for heterogeneity for each of the groups of time-series 
and case-crossover studies only, and both studies combined 
(Table 2). Statistical significance was observed for the 8 time- 
series studies {p< 0.001), but not for the 4 case-crossover stud- 
ies (p=0.65). For the combined studies, we found heterogene- 
ity with the 10 studies for the main analysis (p<0.001) as well 
as all the 12 studies (p<0.001).Thus, we applied random-ef- 



Studies identified from PubMed (n=378), EMBASE (n=441 ), and the Cochran 
Library (n=0) using the keywords "Korea OR Japan OR China OR Mongolia 
OR Thailand OR Malaysia OR Indonesia OR "Vietnam" OR "Papua New 
Guinea" OR Philippines OR Singapore OR Taiwan AND (PM10 OR "air pollu- 
tion" OR particulate) AND mortality" 











186 Excluded duplicates 






633 Studies remaining after duplicates were removed 








585 Excluded according to selection criteria: 

- 68 Studies not presented in English 

- 129 Studies on subjects other than humans 

- 346 Studies with exposures other than to PM10 

- 42 Studies with outcomes other than total or non- 

accidental mortality 






48 Studies reviewed including the whole text 








36 Excluded studies: 

- 8 Studies other than time-series or case-crossover 

analysis 

- 7 Studies with unfit population 

- 3 Studies with insufficient accuracy 

- 18 Studies with overlapping dataset/period/region 






12 Studies finally selected for the analysis 



Figure 1. Selection process of the relevant studies for meta- 
analysis. PM10, particulate matter less than 10 urn in aerody- 
namic diameter. 



fects model to the time-series studies and combined analysis, 
and fixed-effects model to the case-crossover studies. 

Figure 2 shows the forest plots of (A) 1 0 studies for both the 
time-series and case-crossover studies combined and (B) the 
8 time-series and 4 case-crossover studies separately ana- 
lyzed and then the 12 studies all combined. All the results 
showed a significant positive association between PM10 expo- 
sure and non-accidental mortality (RR [95% CI]): (A) 1.0044 
(1 .0033 to 1 .0062); (B) 1 .0057 (1 .0029 to 1 .0086) for the time- 
series studies, 1 .0027 (1 .001 0 to 1 .0043) for the case-crossover 
studies, and 1.0044 (1.0032 to 1.0056) for all the 12 studies 
combined. 

By symmetric Begg funnel plots and non-significant Egger 
tests (p=0.08) on the 1 2 studies, we assessed that this analysis 
was not likely to have publication bias (Figure 3). 
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Table 2. Meta-analysis results for the selected studies 



Study groups 


Study number 


Summary estimate (RR [95% CI]) 


l 2 (%) 


X 2 


X 2 p-value 


Eggers p-value 


Studies from each city 


10 


1.0047(1.0033, 1.0062) 


63.01 


28.2615 


<0.001 


0.15 


Time-series studies only 


8 


1.0057(1.0029, 1.0085) 


92.43 


32.4486 


<0.001 


na 


Case-crossover studies only 


4 


1.0027(1.0010, 1.0043) 


15.09 


1.6214 


0.65 


na 


All selected studies 


12 


1.0044(1.0032, 1.0055) 


54.83 


33.7341 


<0.001 


0.08 



RR, relative risk; CI, confidence interval; na, not applicable. 



City(y) 



Weight (%) 



RR (95% CI) 



1 . Wuhan (2000-2004) 

2. Shanghai (2004-2005) 

4. Inchon (1995-1996) 

5. Hong Kong (1996-2002) 

6. Bangkok (1999-2003) 

7. Guangzhou (2006-2009) 

8. Tianjin (2006-2009) 

10. Anshan (2004-2006) 

11. Seoul (2000-2006) 

12. Kaohsiung (1994-2000) 



15.8582 
17.0926 

4.5379 
12.5563 

5.8723 

7.9028 
16.2830 

3.9797 
15.8496 
H 0.0677 



1.0036(1.0019, 
1.0016(1.0002, 
1.0070(1.0011, 
1.0040(1.0015, 
1.0130(1.0080, 
1.0126(1.0086, 
1.0042(1.0026, 
1.0067(1.0003, 
1.0024(1.0007, 
1.0000 (0.9470, 



1.0053) 
1.0030) 
1.0130) 
1.0065) 
1.0180) 
1.0166) 
1.0058) 
1.0132) 
1.0041) 
1.0560) 



Random-effects model 



100.0000 1.0047(1.0033,1.0062) 



0.9252 1.0015 1.0778 
Observed outcome 



Case-crossover studies 

9. Shanghai (2000-2001) 

10. Anshan (2004-2006) 

11. Seoul (2000-2006) 

12. Kaohsiung (1994-2000) 
Fixed-effects model 



Random-effects model 



City(y) 




Weight (%) 


RR (95% CI) 


Time-series studies 








1. Wuhan (2000-2004) 


■ 


13.9931 


1.0036(1.0019,1.0053) 


2. Shanghai (2004-2005) 


■ 


15.4766 


1.0016(1.0002,1.0030) 


3. Seoul (1999-2004) 


■ 


15.0023 


1.0037(1.0022,1.0052) 


4. Inchon (1995-1996) 




3.2458 


1.0070(1.0011, 1.0130) 


5. Hong Kong (1996-2002) 


■ 


10.3726 


1.0040(1.0015, 1.0065) 


6. Bangkok (1999-2003) 




4.2961 


1.0130(1.0080,1.0180) 


7. Guangzhou (2006-2009) 


H 


5.9898 


1.0126(1.0086,1.0166) 


8. Tianjin (2006-2009) 


■ 


14.4950 


1.0042(1.0026,1.0058) 


Random- effects model 


♦ 




1.0057(1.0029,1.0086) 



0.2806 
2.8202 
13.9830 
I 0.0450 



1.0030(0.9815,1.0250) 
1.0067(1.0003,1.0132) 
1.0030(1.0007,1.0041) 
1.0067(0.9470,1.0560) 
1.0027(1.0010,1.0043) 



1.0044(1.0032,1.0056) 



0.9252 1.0015 1.0778 
Observed outcome 



Figure 2. Forest plot of the studies used in the meta-analysis of particulate matter less than 10 um in aerodynamic diameter expo- 
sure and non-accidental mortality: (A) 1 0 studies for both time-series and case-crossover studies combined; (B) The 8 time-series 
and 4 case-crossover studies, separated and combined. RR, relative risk; CI, confidence interval. 



0.000 



0.041 



0.082 
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0.40 



Figure 3. Funnel plot of all 1 2 of the studies in the meta- 
analysis. 

DISCUSSION 

The estimated association between PMio exposure and non- 
accidental mortality in Asian countries shows a significant pos- 



itive relationship, by time-series studies, case-crossover stud- 
ies, and also by combined analysis. RR ranged from 1.0019 
(1.0005 to 1.0033) to 1.0049 (1.0020 to 1.0077), implying that 
the percentage change of risk per 10 ug/m 3 increase in PMio 
ranged from 0.1 9% to 0.49%. 

Although influenced by the same general trends of modern- 
ization and industrialization, socio-economic changes in Asian 
countries have occurred differently from North American and 
European countries. Asian countries have unique cultural and 
socio-demographic factors as well as different meteorological 
conditions and climate change, and therefore have important 
differences in outdoor air pollution exposure and susceptibility. 

Evaluating the health impact of air pollution in Asian popu- 
lations has been previously attempted by the Public Health 
and Air Pollution in Asia (PAPA) study, where they focused on 
the association between four air pollutants (N0 2 , S0 2 , PMio, 
and 0 3 ) and mortality (non-accidental, cardiovascular, and 
respiratory) in four Asian cities from 1996 to 2004: Bangkok 
(1999 to 2003), Hong Kong (1996 to 2002), Shanghai (2001 to 
2004), and Wuhan (2001 to 2004) [1 3]. Their reported percent- 
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age change or risk per 10 |jg/m 3 increase in PMio, for all four 
cities combined, was 0.55% (0.26 to 0.85), which is larger than 
our estimated effect of 0.33%; however, our study covered 
more regions than the PAPA study and used results from stud- 
ies from a more recent period except for Hong Kong, where 
results reported by the PAPA study were selected. The China 
Air Pollution and Health Effects study on association of PMio 
and daily mortality in 16 Chinese cities which followed the 
PAPA study protocols reported that the pooled estimate for 
the increase in total mortality was 0.35% (0.1 8 to 0.52), which 
is similar to our results [26]. Another recently-published study 
that systemically reviewed the association between air pollu- 
tion and health outcomes in Asia estimated the effect of PMio 
on all-cause mortality from a meta-analysis of eight studies as 
a 0.14% (0.09 to 0.19) change in the fixed-effect summary es- 
timate and 0.27% (0.12 to 0.42) as the random-effect [27]. This 
study is similar to our study with some overlapping studies, 
but it is limited to time-series designs only, while our study 
also includes a case-crossover analysis. We have also used more 
recent studies for Seoul and Bangkok and have included study 
results from the novel regions, Tianjin and Guangzhou. Our 
summary estimate of 0.47% (0.33 to 0.62) is comparable with 
these previous publications in Asia as well as reports from the 
United States (0.29% [0.18 to 0.40] in lag 1, 0.14% [-0.12 to 
0.40] in lag 0 to 1), Canada (0.84% [0.30 to 1.40] in lag 1) and 
Europe (0.33% [0.22 to 0.44] in lag 1, 0.29% [0.14 to 0.45] in 
lag 0to1) [13,28]. 

The World Health Organization has recommended air quali- 
ty guidelines based on extensive scientific review of air pollu- 
tion and its health effects in order to provide useful means to 
protect public health worldwide [29]. As an extension to this 
work, the Thematic Working Group on Air Quality has been set 
up among Asian countries to find practical means for health 
impact assessment of air pollution indigenous to Asian coun- 
tries. Efforts have been made to apply research results to poli- 
cies, and our study may contribute to instituting guidelines for 
PMio exposure levels in Asian countries. 

The toxicological mechanisms for the effect of ambient air 
pollution, including PMio exposure, on humans has been ex- 
plained by pro-inflammatory reaction from increased free rad- 
ical activity by oxidative stress, and increased plasma viscosity 
by stimulation of the mediators for blood coagulation [30-33]. 
Authors have proposed that these systemic reactions may ex- 
tend to non-accidental mortality, including cardiovascular and 
respiratory mortality. Nonetheless, further research is needed 



to elucidate the exact pathological mechanism linking PMio 
exposure to mortality. 

Reports of studies comparing results from both time-series 
and case-crossover analysis have been inconsistent. While sim- 
ilar results have been produced in studies on total suspended 
particles and total mortality [34,35], a study on the association 
between PMio exposure and cardiovascular mortality showed 
that results differed greatly with changing strata lengths [36]. 
Another study on the association between PMio and hospital 
admissions for cardiovascular diseases also showed that results 
were less robust when analyzed by a case-crossover design 
than by using a time-series model [37]. Although with differ- 
ent exposures and outcomes (i.e., diurnal temperature range 
or extreme high temperatures and hospital admissions for car- 
diovascular and respiratory diseases), two other studies that 
compared time-series and case-crossover analyses results re- 
ported that the estimated effects showed wider CIs in the 
case-crossover analyses than the time-series analyses, and the 
superiority of effect sizes was also inconsistent [38,39]. 

In our results, a higher effect size was estimated from the 
time-series studies than the case-crossover studies, and the 
results were statistically significant in both groups. Due to the 
low number of studies used in the meta-analysis, our findings 
may not be the definitive supplement to comparisons of the 
two study designs. However, it is still meaningful to compare 
the two designs, not by examining the superiority of one de- 
sign over another, but by merely comparing the differences in 
the results of the two designs. 

The results from each study varied in population size, mean 
PMio concentrations, and daily mortality count, and therefore 
the size estimate. Low levels of daily PMio levels were reported 
in Hong Kong and Bangkok (51.6 ug/m 3 , 52.1 ug/m 3 ), while 
Wuhan, Shanghai and Anshan in China reported high levels 
(141.8 ug/m 3 , 107.9 ug/m 3 , 1 10.9 ug/m 3 ) which were still be- 
low the US Environmental Protection Agency air quality stan- 
dards (24-h average of 150 ug/m 3 ). The percentage change in 
non-accidental mortality per 1 0 ug/m 3 increase in PMio was as 
high as 1.3% and 1.26% in Bangkok and Guangzhou, and as 
low as 0.16% and even 0% in Shanghai and Kaohsiung. Al- 
though we adjusted for the longitude and latitude in our test 
for heterogeneity and publication bias, it should be noted that 
the differences in each region or city should be taken into con- 
sideration in the interpretation of the meta-analysis results. 

The exposure and outcome measures used in all the studies 
are in agreement, i.e., they all used daily monitored data from 
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designated monitoring stations as PM10 measurement and 
used death data recorded by district institutions based on ICD- 
9 (001 to 799) or ICD-10 (A00 to R99) classification. In terms of 
exposure measurement, use of environmental monitoring data 
may not fully represent individual exposure. The results may 
not be exempt from non-directional misclassification, and it is 
likely that the magnitude of the association may have been 
underestimated. Also, as some studies have found that indoor 
and outdoor PM10 levels are highly correlated [13], outdoor 
measurements, to some extent, can also represent more per- 
sonalized indoor exposure. In terms of outcome measurement, 
studies selected for meta-analyses used either ICD-9 or ICD-10, 
which may have contributed to measurement bias; however, a 
study with a similar design and on a similar subject population 
investigated the influence of change in ICD coding in associa- 
tion between air pollution and mortality, and found that there 
was no significant change in the estimated effects [1 4]. 

However, there are limitations to consider in our study re- 
sults. First, we excluded results with adjustments for co-pollut- 
ants from the selected studies, due to a lack of, or different 
combinations of co-pollutants. As the combined effect of pol- 
lutants may be different in magnitude from the effect of PM 10 
alone and we cannot ignore possible interaction, further in- 
vestigation is recommended with more research on co-pollut- 
ant effects over time. Also, we did not analyze for cardiovascu- 
lar and respiratory mortality in our analysis, although many re- 
searchers have focused on the two outcomes. Compared to 
total or non-accidental mortality, the death count numbers for 
cardiovascular or respiratory mortality are small, and the re- 
sults have been inconsistent in both direction and significance 
on many occasions. On summation of the results in the future, 
meta-analysis of the outcomes is therefore needed. 

The subjects in our meta-analysis were the general popula- 
tion living in the areas studied with ungrouped age range and 
socioeconomic status. However, the effect of air pollution on 
mortality is more conspicuous in susceptible populations such 
as infants, the elderly, and socially deprived groups [19,40]. 
While our results may be generalizable to Asian countries and 
especially to residents of developing Asian cities, applying our 
results to susceptible populations may underestimate the ef- 
fects. 

In conclusion, we found a significant positive association 
between PM10 exposure and non-accidental mortality among 
Asian populations, and this association was confirmed in me- 
ta-analyses of time-series studies, case-crossover studies, and 



both study designs combined. We encourage further investi- 
gations to contribute to the health impact assessment and 
public health management of air pollution in Asian countries. 
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